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Abstract. Within the framework of the well-known curvature models, a fluid lipid 
bilayer membrane is regarded as a surface embedded in the three-dimensional Euclidean 
space whose equilibrium shapes are described in terms of its mean and Gaussian 
curvatures by the so-called membrane shape equation. In the present paper, all 
solutions to this equation determining cylindrical membrane shapes are found and 
presented, together with the expressions for the corresponding position vectors, in 
explicit analytic form. The necessary and sufficient conditions for such a surface to 
be closed are derived and several sufficient conditions for its directrix to be simple or 
CNI ■ self-intersecting are given. 
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O ' 1- Introduction 

OO ' 

By a fluid membrane in this paper we assume a membrane formed in aqueous solution 
by a bilayer of lipid molecules, which are in a fluid state, i.e., the molecules can move 
freely within the monolayer they belong to. The structure of the bilayer is such that 
the hydrophobic tails of the lipid molecules situated in different monolayers face one 
another to form a semipermeable core, while their hydrophilic heads face the aqueous 
solutions on either side of the membrane. It is known that the lipid bilayer is the main 
structural component of all biological membranes, the closed lipid bilayer membranes 
(vesicles) being thought of as the simplest model systems for studying basic physical 
properties of the more complex biological cells. 

The foundation of the current theoretical understanding of the vesicle shapes (see, 
e.g. [20, 30, 31]) can be traced more than thirty years back to the works by Canham 
[7] and Helfrich [14], in which the so-called curvature models have been introduced. 
In these models, the vesicle's membrane is regarded as a two-dimensional surface S 
embedded in the three-dimensional Euclidean space 1R 3 and assumed to exhibit purely 
elastic behaviour described by its mean H and Gaussian K curvatures and two material 
constants associated with the bending rigidity of the membrane. 
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In the model proposed by Helfrich [14], currently referred to as the spontaneous- 
curvature model, the equilibrium shapes of the vesicles are determined by the extremals 
of the curvature (shape) energy 



under the constraints of fixed enclosed volume V and total membrane area A. This 
scheme yields the functional 



Here k c and ka are real constants representing the bending and Gaussian rigidity of 
the membrane, li is the spontaneous curvature (a constant introduced by Helfrich to 
reflect the asymmetry of the membrane or its environment), p and A are the Lagrange 
multipliers (another two real constants) corresponding to the constraints of enclosed 
volume and total area, respectively, whose physical meaning is as follows: p represents 
the pressure difference between the outer and inner sides of the membrane, while A is 
interpreted as a tensile stress or a chemical potential (see, e.g., [5]). The Euler-Lagrange 
equation corresponding to the functional (1) reads 



where A is the Laplace-Beltrami operator on the surface S. This equation, derived in 
[28], is often referred to as the membrane shape equation. It is worth noting that the 
second term in the curvature energy T c does not effect equation (2) since its contribution 
to the overall Lagrangian density is a total divergence as follows from the Liouville's 
form of Gauss's Theorema Egregium (see, e.g., [10]). Actually, for closed membranes 
without edges the integral over the Gaussian curvature K is a topological invariant 
by virtue of the Gauss-Bonnet theorem and therefore it may be disregarded until the 
topology of the membrane remains unchanged. This term, however, plays an important 
role when topological phase transitions are considered (see, e.g., [3, 24]) as well as in 
the theory of fluid membranes with free edges (cf. [9, 36, 37]). 

Later on, another two curvature models have been developed. The first of them is 
the so-called bilayer-couple model suggested by Svetina and Zeks in [33] on the ground 
of the bilayer-couple hypothesis [32] and the related work [34] . The second one is known 
as the area-difference-elasticity model [4, 22, 42]. For the purposes of the present paper, 
however, it is important to underline that all the curvature models mentioned above 
lead to the same set of stationary shapes, determined locally by equation (2), since they 
differ only in global energy terms (see [20, 22, 33]). Of course, the constants involved in 
this equation have different meaning in different models. 

For more than three decades, the study of the equilibrium shapes of the vesicles has 
attracted much attention nevertheless to the best of our knowledge only a few analytic 
solutions to equation (2) have been reported. These are the solutions determining: 
spheres and circular cylinders [28], Clifford tori [29] and toroidal shapes (cf. [20], 
Ch. 8, Sect. 5), circular biconcave discoids [25], Delaunay surfaces [23, 26], nodoidlike 





(1) 



2k c AH + k c (2H + li) (2H 2 -tH - 2K) - 2\H + p = 0, 



(2) 
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and unduloidlike shapes [26], several types of surfaces with constant squared mean 
curvature density and Willmore surfaces (cf. [16, 38, 43]) as well as cylindrical surfaces 
corresponding top ^ [39, 40] andp = 0. The latter case is exceptional since it coincides 
with the prominent Euler's elastica whose typical equilibrium shapes are known for a 
long time [21]. This problem has been studied for more than three centuries in various 
contexts, see, e.g., the recent papers [6, 15, 18, 19, 45]. 

The main goal of the present paper is to present in analytic form all solutions to 
the membrane shape equation (2) determining cylindrical surfaces as well as to provide 
explicit expressions for the corresponding position vectors. In a sense, it might be 
thought of as a completion, from analytic point of view, of the works by Arreaga et al. 
[2] and Capovilla et al. [8] where a purely geometric construction of the equilibrium 
shapes of closed planar loops subject to the constraints of fixed length and enclosed area 
is presented. 

2. Cylindrical equilibrium shapes 

For cylindrical surfaces in H 3 whose directrices are plane curves V of curvature n(s) 
parametrized by their arclength s, the corresponding generatrices being perpendicular 
to the plane the directrices T lie in (see figure 1), the general shape equation (2) simplifies 
and reduces to the single ordinary differential equation 



Indeed, using the standard formulas from the textbooks on classical differential geometry 
(see, e.g., [10, 27]) one can easily find that for a cylindrical surface parametrized in the 
above way H = (l/2)/«(s) and AH = (l/2)d 2 fv(s)/ds 2 in addition to the relation K = 0. 
Substituting the latter expressions in equation (2) one gets immediately to equation (3). 




(3) 



where 





2p 



z 




X 



Figure 1. A slice of the infinite generalized cylinder (left) and its intersection with 
the plane Y = (right). Here, t(s), tp(s) and 9(s) are the tangent vector, slope angle 
and the angle between the position vectors x(0) and x(s), respectively. 
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In what follows, we are interested in real-valued solutions k(s) 7^ const of equation 
(3) possessing smooth derivatives. Once such a solution is known, it is possible to recover 
the embedding x (s) = (x(s),z(s)) E R 2 of the corresponding directrix T in the plane 
1R 2 (up to a rigid motion) in the standard manner. First, recall that the unit tangent 
t (s) = (dx (s) /ds,dz (s) /ds) and normal n(s) = (— dz (s) /ds, dx (s) /ds) vectors to 
the curve T are related to the curvature k(s) through Frenet-Serret formulas [10, 27] 

^ = «(.)„(.), ^ = - K (s)t(s). (4) 
Consequently, in terms of the slope angle cp (s) of the curve T (see figure 1), one has 

«M = ^. ^ = «»(*>«>, ^=sm(„ W ) (5) 
and hence, the parametric equations of the curve T can be expressed by quadratures 

«(.) = / cos („(.)) <U, Z ( S )=/sin (l ,( S ))d S , (6) 

where 

«(a)da. (7) 



Thus, the first problem to solve on the way to determining the cylindrical equilibrium 
shapes of the fluid membranes is to find the solutions of equation (3) in analytic form. 

Fortunately, equation (3) is integrable by quadrature since it falls into the class of 
equations describing conservative systems with one degree of freedom [1]. Indeed, (3) 
can be regarded as the equation of motion of a fictitious particle of unit mass whose 
kinetic, potential and total energies are 

T=-(^)\ U{k) = -k a --hk 2 --(tk, S = T + U, 
2 \ds J ' K ' 8 4 P 2 ' 

respectively. In this setting, k is interpreted as the displacement of the particle while s 
plays the role of the time. The total energy S of this system is conserved and hence 

^) 2 = P(*0, P^) = 2E-\^+ 1 -^ + c-k (8) 

holds on each continuous solution of equation (3), E being the value of its total energy. 
Therefore, the solution of equation (3) can be reduced to the quadrature 

dK = [ dK (9) 

y/2(E - U{k)) J ^/P{k) 

up to a shift of the independent variable s and change of its sign to the opposite one. 
Note that equation (3) and its first integral (8) are invariant under the aforementioned 
transformations of the variable s. Note also that for each solution k = k(s) of equation 
(3) relation (8) implies the existence of a certain value of the variable s at which 
dn/ds = (this matter will be clarified in detail at the beginning of the next Section). 
Without loss of generality, this value may be chosen to be zero due to the translational 
invariance of equation (3). So, hereafter we will always chose dn/ds = at s = 0. 
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Moreover, the specific differential structure of equation (3) allows the integration 
in expressions (6) to be avoided when o ^ 0. Below, one can find a simple alternative 
derivation of this remarkable integrability property, first established in [2] (see also 
[8, 12]). Actually, a direct computation shows that the following identity holds 



2 / dn(s) \ d / d/c(s) , . , 2 , 



■(« 2 (s) - n) [ + K (s)t(s) J = — ( 2-^t(s) - (k 2 (s) - n)n(s) - <rx(s) 



(10) 



and hence, taking into account Frenet-Serret formulas (4), one can represent the position 
vector x (s) of a plane curve of curvature k (s) in the form 

x (s) = -^t (s) - -(k%s) - /i)n (s) + C, 
a as a 

where C is a constant vector, provided that k (s) is a solution of equation (3) with 
a 7^ 0. Then, translating the origin so that x • t = and x • n = — 1/a (k 2 — fi) when 
dft/ds = 0, which is always possible, one gets C = and obtains, taking into account 
the definitions of the tangent and normal vectors as well as the second and the third of 
the relations (5), the following expressions for the components of the position vector in 
therms of the curvature k(s) and slope angle (p(s) 

x(s) = — - ^ cosip(s) H — (^ 2 (s) — /•*) sin</?(s), 

a ds a ( U ) 

z(s) = s smip(s) (k 2 (s) — n) cos ip(s). 

a ds a 

Note, however, that the slope angle (p(s) still remains determined implicitly, via an 
integration (7) of the curvature k(s), which, so far, cannot be accomplished. Note also 
that formulas (8) and (11) lead to the remarkable relation 

= SE + f + 4kW 
a 2 a 



for the magnitude r(s) = \Jx 2 (s) + z 2 (s) of the position vector x(s) found in [2, 8]. 
For <t = 0, the situation is quite similar. In this case, identity (10) implies 

^t( S )-i( K 2 ( S )-/i)n( S ) = C , (13) 

where Co is a constant vector, provided that k (s) is a solution of equation (3) with a — 
and relations (4) hold. Now, choosing (p = dn/ds = at s = and taking into account 
relation (8), the definitions of the tangent and normal vectors as well as the second and 
the third of the relations (5), one can first see that C = (0, — (/t 2 (0) — /i)/2) ^ and 
then rewrite equality (13) in the form 

. k 2 (s) - fi . 2 d/c(s) 

cos " (s) = ^(oP7;' " y( ^" ^ d. • (14) 

Consequently, formulas (6) and (14) imply 

x ( s ) = 1 [K 2 (s)ds- , z(s) = - 9 2/t ! S) ■ (15) 

v; k 2 (0)-hJ KJ /t 2 (0)-Ai' v; /t 2 (0)-Ai K J 
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Again, one integration remains to be done. This time, the square of the curvature n(s) 
has to be integrated. 

Before proceeding with the derivation of the solutions of equation (3) using the 
quadrature (9), it should be remarked that this equation have been regarded in a number 
of papers, see [11, 35, 41, 44], that have not been mentioned yet since they do not concern 
directly the problem considered here. Closest to the subject of the present paper are 
[2, 8] where equation (3) is introduced with the aim to study the equilibria of an elastic 
loop in the plane subject to the constraints of fixed length and enclosed area. In the 
three dimensional case considered here, each such loop determines a directrix T of a 
cylindrical surface whose mean curvature satisfies the membrane shape equation. In 
the foregoing two papers, the authors have succeeded in obtaining a purely geometric 
construction for determination of the curvature of the loop passing through a given 
point of the plane without using explicit expressions for the solutions of equation (3). 
Nevertheless, in our opinion, the knowledge of the solutions to this equation in analytic 
form is of considerable interest for the further exploration of the cylindrical equilibrium 
shapes of the fluid membranes. 

3. Explicit analytic solutions 

First, it should be remarked that the quadrature (9) can be easily expressed in terms of 
elliptic integrals [13] or elementary functions by means of the roots of the polynomial 
P{k) provided that the following observations are taken into account. 

Given a solution k(s) of an equation of form (3) with coefficients /x and a, let E be 
the value of its total energy. Then, bearing in mind that /i, a and E are real numbers, 
it is clear that the corresponding polynomial P(k) has at least two different real roots, 
otherwise the function P(k,(s)) could not take different nonnegative values, as required 
by relation (8) and the assumptions concerning the type of solutions considered, since 
the coefficient at the highest power /t 4 of this polynomial is negative. Thus, in general, 
only two alternative possibilities have to be considered, namely: (I) the polynomial 
P(k) has two simple real roots a, f3 G M, a < ft, and a pair of complex conjugate roots 
7, 5 G C, 5 = 7; (II) the polynomial P(k) has four simple real roots «</?<7<i5eM. 
In the first case, the polynomial P(n) is nonnegative in the interval a < k < (3, while 
in the second one, it is nonnegative in the intervals a < k < (5 and 7 < k < 8. 

It should be noted also that the roots a, (3, 7 and 5 of the polynomial P(k) can be 
expressed explicitly through its coefficients /1, a and E and vice versa. Indeed, after 
some standard algebraic manipulations (cf. [17], Sec. 1.8), one can find the following 
expressions for the roots of the polynomial P(k) 
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where 



uj = 



/i + {3 (3V 2 + Jx) - ii + 2 3 3 2 E) 



- 2 3 3E 



3^3 (3V 2 + v ^)- / u(/i 2 + 233 2 E) 



X = 3 \2 3 E 



(fi 2 + 8E) 2 - 3 2 2/2o- 2 



- a 2 (2^ 



3 3 a 2 ) 



Then, one can denote properly each of the above expressions in accordance with the 
notation introduced in the cases (I) and (II), respectively. Simultaneously, by Vieta's 
formulas one obtains 



due to the absence of a term with k, 3 in the polynomial P(k), and consequently 
H = - (a 2 + j3 2 + 7 2 + a(3 + «7 + ^7) , 
o = - \ (a + (3) (a + 7) {(3 + 7) , 



E 



a(3^ (a + f3 + 7) . 



(16) 

(17) 
(18) 
(19) 



Now, we are in a position to express the arclength as a function of the curvature 
representing the quadrature (9) via elliptic integrals or elementary functions in each 
of the particular cases (I) and (II). Instead of this, however, taking the corresponding 
inverse functions we prefer to give directly the curvature as a function of the arclength 
in terms of the roots of the polynomial P(k). Moreover, solving the integral (7), we 
give explicit formulas for the corresponding slope angles as well. The explicit analytic 
expressions for the solutions of equation (3) and the corresponding slope angles can be 
found in Theorem 1, for the case (I), and in Theorem 2, for the case (II). 

Theorem 1 Given /i, a and E, let the roots a, f3, 7 and 5 of the respective polynomial 
P(k) are as in the case (I), that is a < j3 G R, 7, 5 G C, 5 = 7, and let 77 = (7 — 7)/2i. 
Then, except for the cases in which (3a + (3)(a + 3/5) = r\ = 0, the function 

(Af3 + Ba) - (A(3 - Ba) cn(ws, k) 



Kl (S) (A + B) - (A - B) cn( MS , k) 
of the real variable s, where 

A= yV + (3a + /3) 2 , 



B 



^47? 2 + (a + 3/3) 2 , u = j^/AB , 



y/2' 



4r] 2 + (3a + /3)(a + 3(3) 



1- 1 V /V ' =, 

\ ^[Ar] 2 + (3a + (3) (a + 3f3)] 2 + lQr] 2 ((3 - a) 2 



(20) 

(21) 
(22) 



takes real values for each s G R and satisfies equation (3). This function is periodic if 
and only ifr/^0 or r\ = and (3a + (3)(a + 3(3) > 0, its least period is Ti = (4/w) K (k), 
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and for a 7^ its indefinite integral tpi (s) such that (pi (0) = is 

<Pi{8)- A _ B s+ 2u{A _ B) HI 4AB ,&m(us,k),k 



<\ - I I /, , (A - B) 2 m(us,k) X 
H . = arctan \ k z J \ — — — -, 

W*' + V V ^ ^"'^ 

In £/ie cases m which (3a + /3)(a + 3/3) = n = 0, £ne function 

where C = ot when 3a + /3 = and C = /? w/ien a + 3/3 = 0, satisfies equation (3) for 
each s6l and its indefinite integral ip 2 (s) such that ip 2 (0) = reads 

V^2 (s) = (s — 4 arctan (£s) . (25) 

Proof. Let us begin with the simpler case (3a + (3) {a + 3/3) = n = 0. Here, on account 
of the definition of r] and relations (16)— (19), a straightforward computation shows that 
the function (24) is the derivative of the function (25) and satisfies equation (3). The 
relation <p 2 (0) = is obvious. 

Next, let us exclude the cases in which (3a + (3) {a + 3(3) — rj — 0. Now, the 
condition a < (3 G R, the definition of 77 and expressions (21) and (22) imply also that 
T],A,B,u,k G R, AB 0, u > and < fc < 1. Hence, the function (20) is real- 
valued when s G R. Substituting the function (20) into equation (3) and taking into 
account relations (16)-(19), one can easily verify that the latter equation is satisfied. 
When rj 7^ 0, expression (22) implies < k < 1 and therefore the function (20) 
is periodic because the function cn (us, k) is periodic, with least period T\. When 
77 = but (3a + (3)(a + 3(3) 7^ two alternative cases are to be considered, namely 
(3a + (3)(a + 3(3) > and (3a + (3)(a+ 3(3) < 0. In the first case, expression (22) leads 
to k — 0, which means that cn(ws, k) = cos(ws) and hence the function (20) is periodic 
again. However, if (3a + (3)(a + 3(3) < 0, then expression (22) implies k — 1 meaning 
that cn(us,k) = sech(ns), and hence the function (20) is not periodic. Thus, having 
considered all possible cases, we may conclude that the function (20) is periodic if and 
only if rj 7^ or rj = and (3a + (3) (at + 3(3) > 0. Finally, differentiation of expressions 
(23) with respect to the variable s yields (20). The relation tpi (0) = is obvious. ■ 
Theorem 2 Given fi, a and E, let the roots a, (3, 7 and 5 of the respective polynomial 
P(k) are as in the case (II), that zsa</3<7<5GlL Consider the functions 
, n c (5 -a) (5 -(3) 

^^ = 6 ~ (*-/?) + (/?- a) sn' («*,*)' ( } 

k <*\ <? ■ (7-/3)0*-/?) (97) 
K 4 (s) = (3 + — r — — — , (27) 

(0 — (3) — (0 — 7) sn 2 (us, k) 
of the real variable s, where 
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Then, both functions (26) and (27) take real values for each sGffi and satisfy equation 
(3), they are periodic with least period T 2 = (2/u)K(k) and their indefinite integrals 
ip 3 (s) and ip4 (s), respectively, such that ip 3 (0) = ip± (0) = are 

(p 3 (s) = 5s — II I - _ - , am(ws, k), k \ , (29) 

V?4 (s) = I3s — — — -IT ( - — ^, ain(«s, k),k) . (30) 
u \o — p ) 

Proof. It is easy to see that the condition «</?<7<ielR and expressions (28) 
imply also that u, k G R, u > and < /c < 1. Therefore, both functions (26) and 
(27) are real-valued when s e R. Substituting each of the above functions into equation 
(3) and taking into account relations (16)— (19) , one can easily verify that they satisfy 
it. Evidently these functions are periodic due to the fact that the function sn 2 (us, k) 
is periodic, with least period T 2 = (2/u)K(k), since u > and < k < 1. Finally, 
differentiation of expressions (29) and (30) with respect to the variable s yields (26) and 
(27), respectively. The relations (f 3 (0) = y?4 (0) = are obvious. ■ 

Suppose now that a = 0. Under this assumption, in the case (I), formulas (16) and 
(18) imply (3 = —a > 0. Then, according to Theorem 1, B = A = 2^rf + a 2 and so 



1 / a 2 

Ki(s) — acn(us,k), u=-^/r] 2 + a 2 , k = \ ^— — -, (31) 

V ^ 

cf. formulas (20) - (22). Consequently 

K 2 1 (s)ds = 2 V 'a 2 + r] 2 E (am(ws, k), k) - rfs. (32) 



Note also that formulas (31) and (17) imply 

«i(0) = a, fi = ^(a 2 -r ] 2 ). (33) 

In the case (II), the assumption a = and formulas (16) and (18) imply 5 = —a > 0, 
7 = —f3 > 0. Hence, according to Theorem 2, cf. formulas (26)-(28), we have 

and 

1 — /csn 2 (us, k) . 1 + k sn 2 (us, k) 

1 + k sn 2 (us, k) ' ' 1 — k sn 2 (-us, k) 

Then, using Gauss's transformation u — u (1 + k), k = 2\fkj (k + 1), we obtain 

«3(s) = a dn(-us, k), u — —-a, k = \/ a 2 — /5 2 , (34) 

/t 4 (s) =7 — — = -adn(us + K(k),k) = -k 3 (s + m _1 K(A;))- (35) 

dn(ws, A;) 

Now, observing relation (35), we arrive at the conclusion that the Euler's elastic curves of 
curvatures k 3 (s) and k 4 (s) determined by formulas (34) and (35), respectively, coincide 
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up to a rigid motion in the plane M?. Therefore, to complete the present task it suffices 
to consider only the curves of curvature given by formula (34). In this case, we have 

J k 3 (s) ds = — 2aE(am(-us, k), k) (36) 

and 

« 3 (0)=a, /i=i(a 2 + /3 2 ). (37) 

due to the relations (34) and (17). 

Thus, having obtained in explicit form the solutions of equation (3), i.e., the 
curvatures, the corresponding slope angles and the integrals of the squared curvatures 
(in the cases in which a = 0), we have completely determined in analytic form the 
corresponding directrices T (up to a rigid motion in the plane R 2 ) through the parametric 
equations (11) when a ^ or (15), with the supplementary relations (31)-(33) or (34), 
(36) and (37), when a — 0. Several examples of directrices of cylindrical equilibrium 
shapes corresponding to solutions to equation (3) of form (20), (26) or (27) are presented 
in figure 2 and figure 3 for a ^ and in figure 4 for a = 0. 




Figure 2. Directrices of some closed cylindrical equilibrium shapes whose curvatures 
are solutions to equation (3) of form (20) with u = 1 and: a) /z = 1.908, E — 0.146; 
b) n = 0, E = 0.211; c) jj, = 1/3, E = 0.407; d) /z = 1/3, E = 0.563. 




Figure 3. Directrices of some closed cylindrical equilibrium shapes whose curvatures 
are solutions to equation (3) of forms (26) (top) and (27) (bottom) with [i = 3, a = 1 
and: a) E = 0.0850056; b) E = 0.0849733; c) E = 0.0850046. 
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4. Closure conditions 



Hereafter we are interested in directrices T that close up smoothly meaning that there 
exists a value L of the arclength s such that x(0) = x(L) and t(0) = t(L). The later 
property of such a smooth closed directrix T and the definition of the tangent vector 
imply that there exists an integer m such that y?(X) = <£>(0) + Irrrn where (p(s) is the 
corresponding slope angle. Since throughout this paper we always choose (p(0) = this 
means that the condition t(0) = t(L) implies (f(L) = 2mir. 
Under the above assumptions, expressions (11) show that 



x(0) 



2 dm 



a 



ds 



a 



iK\0)-n)\, x(L) 



2 dn(s) 
a 



ds 



a 



: {k\L)-^) 



Consequently, the equality x(0) = x(L) yields 



dn (s) 



ds 



ds 



k(L) = ±k(0), 



s=0 



which, on account of relation (8), implies that L is a period of the curvature n(s), that 
is L = nT where n is a positive integer and T is the least period of the function k(s). 
Since (p(nT) = rup(T), as follows by formula (7), then 2m7r = <p(L) = tp{nT) = mp(T) 
and hence 

/mN 2miT , . 

<p(T) = (38) 

n 

Thus, in the cases when o ^ 0, relation (38) is found to be a necessary condition for a 
directrix V to close up smoothly. Apparently, it is a sufficient condition as well. 
Straightforward computations lead to the following explicit expressions 



yiy l! u(A-B) y ' u(A-B) \ AAB 



<p 3 (T 2 ) = -K (k) + 2 n — £, k 

u u \ o — p 

ip, (T 2 ) = (k) + 2^n ( p^, k 
u u \p — 



k 
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= 0, «(L) = k(0), x (L) - x (0) = 

=L 



for the angles of form (23), (29) and (30), respectively. These expressions and the closure 
condition (38) allow to determine whether a curve of curvature (20), (26) or (27) closes 
up smoothly or not. 

In the cases when a = 0, relations (14) and (15) show that 
dn (s) 
ds 

are the necessary and sufficient conditions for the respective directrices T to close up 
smoothly. Let us recall that relations (14)-(15) were derived assuming (p = dn/ds = 
at s = 0. The aforementioned conditions mean first that L = nT where n is a positive 
integer and T is the least period of the function k(s) and, consequently, that 

fj,T= [ n 2 (s)ds, (39) 
Jo 

in view of the first of relations (15) and due to the fact that k(s) is a periodic function. 
Now, using formulas (32) and (36), one can easily see that 

/ 1 k j (s) ds = 8 ^a 2 + r] 2 E (k) - tfT u [ * k\ (s) ds = -4aE(jfe) 
Jo Jo 

and then, expressing the coefficient fi from formulas (33) or (37), to rewrite the closure 
condition (39) for the curves of curvatures given by formulas (31) or (34) in the form 

2E(k) - K(k) = 

or 

2E(jfc) - (2 - P)K(jfe) = 0, 

respectively. 

An interesting property of the curves of curvatures k 3 (s) and k^s) is observed in 
the cases in which a ^ (see figure 3). Namely, if one of these curves closes up smoothly, 
then so does the other one. To prove this let us first note that the solutions k 3 (s) and 
«4(s) correspond to the case (II) when the polynomial P(n) has four real roots. Without 
loss of generality, these roots can be written in the form 

a = —3q — v — 2w, j3 = q — v — 2w, 7 = q — v + 2w, 5 = q + 3v + 2w, (40) 

where q, v and w are three arbitrary positive real numbers. The main advantage of this 
parametrization is that it preserves the order of the roots of the polynomial P(k), i.e., 
a < (3 < 7 < 5 for any choice of the parameters q, v and w, which allows one to deal 
freely with them. Using these parameters, it is easy to find that 

meaning that ip = const. This constant can be determined by evaluating the function 
tp for any values of the parameters q, v and w, say q = v = w, which gives 

^ 4 (T 2 ) - y? 3 (T 2 ) = 4tt. (41) 

This relation and the closure condition (38) do imply that the foregoing two curves close 
up simultaneously. 
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5. Self-intersection 

In what concerns the vesicle shapes, of special interest are solutions to equation (3) 
with a 7^ that give rise to closed non-self-intersecting (simple) curves. A sufficient 
condition for such a closed curve to be simple is /i < 0, which is discussed in [8]. It is 
also mentioned therein that the closed curves satisfying condition (38) with m ^ il or 
n — 1 are necessarily self-intersecting. In this Section, the case /i > is considered and 
several new sufficient conditions are presented for a closed curve of the foregoing type 
meeting a closure condition of form (38) with m — ±1 and n > 2 to be simple or not. 

It is convenient to treat the problem of self-intersecting in terms of the magnitude 
r(s) of the position vector x(s) and the angle 9(s) between the position vectors x(0) 
and x(s). Assuming that the angle 9(s) is positive when measured counterclockwise 
from the vector x(0) to the vector x(s) and negative otherwise and taking into account 
expressions (11) we obtain the relations 

x (s) = — sgn(^ (0))r (s) sin 9 (s) , z (s) = sgn(z (0))r (s) cos 9 (s) , (42) 

jgf) = " 2 (s)-J' . (43 ) 
as ar 2 (s) 

The following observation is crucial for the rest of the present study. 

Lemma 1 Let V be a smooth closed curve whose curvature k(s) is a solution of equation 
(3) with a 7^ of form (20), (26) or (27). Let T be the least period of the function k(s) 
and let the corresponding slope angle meets a closure condition of form (38) with m — ±1 
and n > 2, i.e., ip (T) = ±2ir/n. Then, the curve T is self-intersecting if and only if 
there exists s G (0,T/2) such that 9(s ) = 9(0) or 9(s ) = 9(T/2). 

Proof. First, let us recall that each such curve has n axis of symmetry (cf. [2, 8]), any 
position vector x(«T/2), i = 0, 1, . . . , 2n being along one of them. Consequently, any 
point of the curve whose position vector is collinear with one of the foregoing vectors 
lies on an axis of symmetry of the curve T. Note also that on account of formulas (23), 
(29) or (30) the closer condition may be written in the form cp (T/2) = ±7r/n. It should 
be mentioned as well that formulas (20), (26) and (27) imply dft/ds = at s = iT/2. 

Let the curve T be such that 9(sq) = 9(0) or 9(sq) = 9(T/2) at some Sq G (0,T/2). 
Then, expressions (11), (42) and the closer condition imply 9 (T/2) = ±(n/n + In), 
where / is an integer, meaning that the position vector x(so) is along a certain axis of 
symmetry of the curve T. Therefore this curve self- intersects since it passes through two 
different points lying on one and the same axis of symmetry. 

Next, suppose that the curve T is such that 9(s ) ^ 9(0) and 9(s ) ^ 9(T/2) for 
each s G (0,T/2). Then, the same holds true in the next interval (T/2,T) since the 
curve T is symmetric with respect to the axis corresponding to the angle 9(T/2), that 
is along the vector x(T/2), and so on - up to the last interval (nT — T/2,nT). In this 
way, we arrive at the conclusion that the considered curve does not pass twice through 
any one of its axes of symmetry and therefore it is simple because it is simple in the 
interior of each of the foregoing intervals too as is evident from relation (12). ■ 
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Theorem 3 Let T be a smooth closed curve, which meets the assumptions of Lemma 1. 
Then: 

(i) the curve Y is simple if k 2 (s) — ji ^ for s G [0, T/2]; 

(ii) the curve T is self-intersecting if the equation k 2 (s) — /i — has exactly one 
solution for s G [0, T/2}. 

Proof, (i) The condition k 2 (s) -/i ^ and expression (43) imply d9(s)/ds ^ for 
each s G [0, T/2], meaning that 9(s) is a strictly increasing or decreasing function in 
this interval. Consequently, there does not exist s G (0, T/2) such that 9(s ) = 9(0) or 
0(sq) = 9(T/2), and hence, by virtue of Lemma 1, the corresponding curve T is simple. 

(ii) Let s M be the only value in [0, T/2] such that k 2 (s^) — /i — 0. First, let 
s M G (0, T/2). Since k(s) is strictly increasing in this interval, the signs of function 
k 2 (s) — /J, in the intervals (0, s M ) and (s^, T/2) are different. Then, expression (43) 
implies that the function #(s) has an extremum at s M . Suppose, #(s) has a maximum 
at s M . Then, when s increases from to T/2, the function 9(s) increases from 9(0) = 
to a certain value 9(s^) = 9 max > 0, and after that, decreases from 9 max to 9(T/2). 
If 9 (T/2) is negative, then there exists si G (s^, T/2) such that 0(si) = 9(0) meaning 
that the curve T is self-intersecting. If 9 (T/2) is positive, then there exists s 2 G (0, s M ) 
such that 0(s2) = 9(T/2) meaning that the curve T is self-intersecting again. The case 
when 9(s) has a minimum at s M is similar. Next, let s M = or s M = T/2. Then, 
expressions (11) imply x(0) = x(T) = or x(T/2) = x(3T/2) = 0, respectively. On 
the other hand, t(0) 7^ t(T) in the first case because <p(T) = 2n/n, while in the second 
case t(T/2) ^ t(3T/2) since ip(T/2) = n/n and <^(3T/2) = 3vr/n. Therefore, the 
corresponding curves T is self-intersecting. ■ 

Corollary 1 Under the assumptions of Lemma 1, the curve V is self-intersecting if its 
curvature k(s) is such that the respective polynomial P(k) has onlyreal roots. 

Proof. In case (I), according to Theorem 1, real roots are achieved only if 77 = 
meaning that 7 = 5 = —(a + (3)/2. Then, formula (17) implies (a 2 — fi)(/3 2 — [/,)< 
provided that (3a + (3)(a + 3/3) > which is the necessary and sufficient conditions 
for periodicity of the curvature k±(s) when r\ = (see Theorem 1). Consequently, the 
equation k\(s) — /j, = has only one solution for s G [0,T/2] since «i(0) = a and 
Ki(T/2) = f3. Then, according to Theorem 3 (ii), the curve is self-intersecting. 

In case (II), the roots can be written in the form (40) and the angles <fs(T) and 
724 (T) may be thought of as functions of the parameters q, v and w. Differentiating 
relation (41) with respect to the parameter w, one obtains 

dy? 3 (T) = dy? 4 (T) = v 2 - q 2 

dw dw w ( q + v + w ) ^(g + w) (v + w) 

The function E(k) > 1 since < k < 1 and therefore the above derivatives are positive 
when v > q, negative when v < q and equal to zero when v = q. If v > q, then ^ps(T) 
and (pi(T) are increasing functions of the variable w and meet the inequalities 

V? 3 (T) < -2tt, <p A (T) < 2tt 
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since 

lim ip 3 (T) = -2?r, lim ip 4 (T) = 2tt. 

The first inequality implies that if the curve corresponding to the angle <f3(T) closes up, 
then it necessarily self-intersects since m < — 1 in the respective closure condition (38). 
Consider now the curve, corresponding to (fi{T). It is easy to see that 

(7 2 - fi){5 2 - //) = 4[{q + vf + 4vw][q 2 - 3v 2 - 2qv - 4vw] < 

since in this case q 2 — 3v 2 < 0. Consequently, the equation k 2 (s) — \x = has only 
one solution for s e [0, T/2] since ^(0) = 7 and k 4 (T/2) = 5. Therefore, according to 
Theorem 3 (ii), the curve corresponding to (p<i(T) is self-intersecting. Thus, in the case 
v > q both curves are self-intersecting. Similar arguments hold in the case v < q and 
lead to the same conclusion. Finally, if v — q, then the functions y?3(T) and y?/i(T) do 
not depend on w. Evaluating them at q — w, one obtains 

V9 3 (T) = -2tt, V4 (T) = 2tt 

and hence the two curves are self-intersecting as well. ■ 

6. Concluding remarks 

In this paper all solutions to the membrane shape equation (2) determining cylindrical 
surfaces are presented in analytic form in Theorems 1 and 2. Explicit analytic 
expressions for the corresponding slope angles in the case a 7^ are also given in 
these Theorems. Explicit expressions for the integrals of the squared curvatures are 
obtained in the case a = 0. In this way, we have completely determined in analytic 
form the corresponding directrices T (up to a rigid motion in the plane M 2 ) through 
the parametric equations (11) or (15) depending on whether a 7^ or a = 0. The 
necessary and sufficient conditions for the foregoing cylindrical surfaces to be closed are 
derived in Section 4. One necessary and sufficient and three sufficient conditions for 
their directrices to be simple or not are presented in Section 5. These conditions show 
that simple directrices could be achieved only in case (I), i.e., when the polynomial P(n) 
has two simple real roots and a pair of complex conjugate roots. 
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